Long-term analysis of the Störmer-Verlet method for Hamiltonian systems with a solution-dependent high frequency

نویسندگان

  • Ernst Hairer
  • Christian Lubich
چکیده

The long-time behaviour of the Störmer–Verlet–leapfrog method is studied when this method is applied to highly oscillatory Hamiltonian systems with a slowly varying, solution-dependent high frequency. Using the technique of modulated Fourier expansions with state-dependent frequencies, which is newly developed here, the following results are proved: The considered Hamiltonian systems have the action as an adiabatic invariant over long times that cover arbitrary negative powers of the small parameter. The Störmer–Verlet method approximately conserves a modified action and a modified total energy over a long time interval that covers a negative integer power of the small parameter. This power depends on the size of the product of the stepsize with the high frequency.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Dynamic Harmonic Analysis of Long Term over Voltages Based on Time Varying Fourier series in Extended Harmonic Domain

Harmonics have become an important issue in modern power systems. The widespread penetration of non-linear loads to emerging power systems has turned power quality analysis into an important operation issue under both steady state and transient conditions. This paper employs an Extended Harmonic Domain (EHD) based framework for dynamic analysis of long term analysis over voltages during the tra...

متن کامل

Second Order Conformal Symplectic Schemes for Damped Hamiltonian Systems

Numerical methods for solving weakly damped Hamiltonian systems are constructed using the popular Störmer-Verlet and implicit midpoint methods. Each method is shown to preserve dissipation of symplecticity and dissipation of angular momentum of an N -body system with pairwise distance dependent interactions. Necessary and sufficient conditions for second order accuracy are derived. Analysis for...

متن کامل

Ergodicity and the Numerical Simulation of Hamiltonian Systems

We discuss the long-time numerical simulation of Hamiltonian systems of ordinary differential equations. Our goal is to explain the ability of symplectic integration schemes such as Störmer-Verlet to compute accurate long-time averages for these systems in the context of molecular dynamics. This paper introduces a weakened version of ergodicity that allows us to study this problem. First, we de...

متن کامل

Exponential Integration for Hamiltonian Monte Carlo

We investigate numerical integration of ordinary differential equations (ODEs) for Hamiltonian Monte Carlo (HMC). High-quality integration is crucial for designing efficient and effective proposals for HMC. While the standard method is leapfrog (Störmer-Verlet) integration, we propose the use of an exponential integrator, which is robust to stiff ODEs with highly-oscillatory components. This os...

متن کامل

Symplectic Time-Stepping for Particle Methods

This paper surveys some of the fundamental properties of symplectic integration schemes for classical mechanics and particle methods in particular. The widely used Störmer-Verlet method is discussed in detail and implications of conservation of symplecticity on long term simulations are outlined. The second part of the paper describes the application of a Lagrangian particle method and the Stör...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Numerische Mathematik

دوره 134  شماره 

صفحات  -

تاریخ انتشار 2016